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The probability distribution is found for the 
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tioned mobile radios in a wireless network 
for two representative deployment scenar- 
ios: (1) the mobile locations are uniformly 
distributed over a rectangular area and 
(2) the x and y coordinates of the mobile 
locations have Gaussian distributions. It 
is shown that the shapes of the link dis- 
tance distributions for these scenarios are 
very similar when the width of the rectan- 
gular area in the first scenario is taken to 
be about three times the standard deviation 
of the location distribution in the second 



scenario. Thus the choice of mobile loca- 
tion distribution is not critical, but can be 
selected for the convenience of other as- 
pects of the analysis or simulation of the 
mobile system. 
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1. Introduction 



The probability that a link between two mobile radios 
has sufficient signal-to-noise ratio for acceptable trans- 
mission quality or reliability is, other factors being 
equal, the probability that the link distance d is less that 
some value R, where R is termed the transmission 
range: 

Pr{Link is good} = ¥r{d < R] = F d (R). (1.1) 

The function F d (-) in Eq. (1.1) is the cumulative proba- 
bility distribution function (cdf) for the link distance. 

Assuming that different links fail independently, the 
quantity F d (R) can be taken as the probability of suc- 
cess (acceptable transmission quality) in a binomial trial 
in which two link endpoints are selected; if the trial is 
repeated N times, then an estimate of the number of 



good links is NF d (R). Also, the probability that multi- 
hop communication paths are reliable can be related to 
the individual link reliabilities. For these and other rea- 
sons, the cdf for the link distances in a mobile radio 
system is an important quantity [1,2]. 

There is an infinite number of potential scenarios in 
which locations are selected for the different mobile 
radios. In this paper, in order to lay the goundwork for 
further analysis of mobile radio systems, a random selec- 
tion of mobile locations is assumed, and the cdf of the 
link distances is found for two simple but fundamental 
scenarios: (1) a rectangular deployment area in which 
mobiles are uniformly distributed and (2) a deployment 
in which the x and y coordinates of the mobile locations 
have Gaussian distributions. 
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2. Uniform Distribution of Link Distances 
in a Rectangular Area 

2.1 Assumptions and Formulation of the 
Derivation 

Let the positions of the mobile users (referred to as 
"mobiles") be distributed randomly in a rectangular area 
with dimensions D\ and Z) 2 , as illustrated in Fig. 1, in 
which we have assumed D { < D 2 without loss of gener- 
ality. The x t and y t coordinates of mobile i have the 
uniform distributions given by the probability density 
functions (pdfs) p x (a) and p y (j3), respectively, where 



The link distance between mobiles i and j is defined 

as 

dij = y/fr-xjf + ty-yj) 2 = V(Ax) 2 + (Ay) 2 (2.2) 

where, as illustrated generically in Fig. 2, the differ- 
ences Ax = Xt — Xj and Ay = y t — jj are independent and 
have the pdfs given by 



D, - \a\ 



PAx(a) = ■ 



D] 



.0, 



, \a\<D } 



otherwise 



(2.3a) 



Px{a)- 



W ,\a\^ Dl 



0, otherwise 



(2.1a) 



P>(P) = - 



D 2 



1/31 : 



\Th 



(2.1b) 



0, otherwise 



We assume that the x and y positions of any two mobiles 
are selected independently. 




and 



PAy(P)- 



'^.I0l*ft 



(2.3b) 



.0, 



otherwise 



and where the absolute values of the differences 
IAxI = \x t — Xj\ and lAyl = \y t — yj\ are independent and 
have the pdfs given by 



r 2(Di ~ a) 



P\A.x\(a) = ■ 



and 



P\Ay\(p) = 



D 



Lo, 



, < a < A 
otherwise 



(2.4a) 



pV^o^z> 2 



(2.4b) 



Lo, 



otherwise 



The cumulative probability distribution function for the 
distance between two mobiles therefore is formulated as 



F,( r ) = Pr{^< r }=Pr{V(x / -X;) 2 + (y / -y,) 2 <y} 
= Pr|V(Ax) 2 + (Ay) 2 <y| = Pr|VlAxl 2 + IAyl 2 <y 

// da d/3/?iAxijAyi(a, )3) 

(a,/3)GA 



(2.5a) 



= 1 - J J da dPp^ Ay] (a, (3) (2.5b) 
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Fig. 2. Pdfs for location, difference, and absolute value of difference. 



where p^^ia , P) = P\Ax\(a)p\A y \(f3) denotes the joint 
pdf of the absolute values of the x and y differences and 
A denotes th e domain of integration, illustrated in Fig. 3, 
such that Va 2 + /3 2 < y while both < a < D { and 
0<j6<Z) 2 , or < a < mintA, y} and 
< j8 < min{Z) 2 , Vy 2 - a 2 }. Using the pdfs of Eqs. 
(2.4a) and (2.4b), Eq. (2.5a) becomes 

/*mm{D\,y) 

F„(y) = / da 






/•min{l,yDi} /rninj 1, ■\ly 2 —D\u 2 iD 2 1 

= 4 / dK(l - h) / dv(l - v) 

(2.6b) 

/•min{l,f} /-min{l,f ^-h 2 ] 

= 4 / dK(l -u) I dv(l - v) 

(2.6c) 

in which we define the normalized variable £ = ylD x 
and the area shape parameter £ = Di/D 2 ^ 1 • The evalu- 
ation of this double integral is facilitated by considering 
different intervals for the value of y. For y < , of 
course, the integral equals zero. For y > VZ) 2 + D\, the 
double integral equals one. Similarly, Eq. (2.5b) be- 
comes 



«-/;-/>5W(-5)('-d 

/ dw(l — u) J 



= 1-4/ dw(l 
with the lower limits 



(2.7a) 
dv(l - v) (2.7b) 



U = 



L\=- 



L 2 = 



L' 2 = 



0, 



< y < D 2 



Vy 2 - D\, D 2 < y < VD 2 + D 



0, 



o<fsr 



vf - r 2 , r'<f^vi + r 



0, < y ^ D, and a > y 



Vy 2 - a 2 , D, < y < VD 2 + D 



0, 0<£<landM>£ 



^vf - w 2 , i<^<vTTr 



(2.7c) 



(2.7d) 
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Fig. 3. Domain of integration for the pdf. 



2.2 Representative Results for the cdf 

In Appendix A, it is shown that the cdf for the link 
distance between two mobiles that are randomly posi- 



tioned in a rectangular area is given by Eq. (2.8). For a 
square area with D t =D 2 = D, or £= 1, the cdf reduces 
to Eq. (2.9). 



F„(y =#),) = 



0, 



et 2 \ket 2 - tea + &+*], 



uvw^ (2f + 1) - \a%e + ete - o 

+ 2tthm\\l£), 



cVF^{2e+D-\c\t+2e-'i) 



0<£<1 



+ 2££ 2 {sin- 1 (l/£) - cos-^l/^)}, c' =s £< viTr 1 

i, VTTr 5 ^ 



(2.8) 



F d (y=£D) = 



0, 



£<0 



f 2 ^ 2 -!f+iT), 0<£<1 



fV?^T(2^+l)-(^ 4 + 2^-|) 

+ 2£ 2 [shr 1 (l/£) - cos-'U/f )], 1 < £< V2 

1, V2<£ 



(2.9) 
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Example plots of Eqs. (2.8) and (2.9) are shown in Fig. 
4. For example, note from Fig. 4 that the median link 
distance (the value of y for which the cdf equals 0.5) is 
approximately y = d me d ~\D f° r me case of £ = L In 
fact, solving F</(£) = 0.5 numerically for £=1 yields 
£med = 0.5120. Additional median values for this distri- 
bution are given in Table 1 for different values of £. 



Pd (y=^D ] ) = - 



2.3 The pdf and Mode for the Link Distance in a 
Rectangular Area 

The probability density function pd(y= £D\) for the 
link distance in a rectangular area is found by differenti- 
ating the cdf in Eq. (2.8) to obtain Eq. (2.10). For the 
special case of D } =D 2 = D or f = 1, Eq. (2.10) be- 
comes Eq. (2.1 1). Example plots of these functions are 
shown in Fig. 5. 




Fig. 4. Plot of the link distance cdf for a rectangular deployment area (D\ = £D 2 ^ D 2 ). 

tf [2tf 2 - 4£(1 + f) + 2ir], 0<£<1 

4#vi 7 =T+4f 2 £V£ 2 -r 2 

- 2£(£ 2 £ 2 + 1 + £ 2 ) 
+ 4tf {siiT'U/fl - cos-'d/tf )}, r 1 =s £< viTr 5 

0, otherwise 



(2.10) 



p rf ( r =£D) = - 



2£(£ 2 -4£+tt), 
8^VF^T-2^ 2 + 2) 



0<£<1 



+ 4£{sirT 1 (l/£) - cos-'U/f)}, 1 < £< V2 

0, otherwise 
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Fig. 5. Plot of the link distance pdf for a rectangular deployment area (D\ = £D 2 ^ D 2 ). 



Table 1. Median values of link distances for a D\ X D 2 rectangular 
area, normalized by D\ < D 2 



Table 2. Mode of the link distances for a D\ X D 2 rectangular area, 
normalized by D\ < D 2 



£ = D,ID 2 



£med = 7mcd/£>l 



£ = D,ID 2 



Smode — 



D 1 



1.00 
0.95 
0.90 
0.85 
0.80 
0.75 
0.70 
0.65 
0.60 
0.55 
0.50 
0.45 
0.40 
0.35 
0.30 
0.25 



0.5120 
0.5254 
0.5401 
0.5563 
0.5743 
0.5943 
0.6170 
0.6428 
0.6725 
0.7072 
0.7486 
0.7990 
0.8625 
0.9465 
1.0666 
1.2453 



From differentiation of the pdf and solving the result- 
ing quadratic equation, the mode of the distribution is 
found to be 



1.00 
0.95 
0.90 
0.85 
0.80 
0.75 
0.70 
0.65 
0.60 
0.55 
0.50 
0.45 
0.40 
0.35 
0.30 
0.25 



0.4786 
0.4908 
0.5034 
0.5165 
0.5299 
0.5439 
0.5582 
0.5730 
0.5882 
0.6037 
0.6196 
0.6357 
0.6521 
0.6687 
0.6855 
0.7023 



smaller than the median values in Table 1, indicating a 
significant amount of skew in the distribution, which 
can be observed in the pdf plots in Fig. 5. 



> ymode 2(1 + / 4(1 + D 2 77 



Example values of the mode for different values of £ are 
given in Table 2. The mode values in Table 2 are 
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3. Distribution of Link Distances for 
Gaussian-Distributed Coordinates 

3.1 Derivation of the Link Distance pdf and cdf for 
Gaussian-Distributed Locations 

Instead of assuming that the mobiles are randomly 
located in a rectangular area, we now assume that the x 
and y coordinates of the mobile locations have Gaussian 
distributions. That is, we assume that the pdfs of the x 
and v coordinates are independent and have the follow- 
ing pdfs: 



p x (a) = } r — e"" 2 ' 2 ^, - oo < a < oo (3.1a) 



ity, we assume that v A = A<x 2 , where A is an area shape 
parameter, with A < 1 . The joint pdf of the coordinates 
is given by 



^•»-sss-*HteH£) 



(3.2a) 



Note that the joint pdf in Eq. (3.2a) is the special case 
of the bivariate Gaussian pdf with uncorrected random 
variables (RVs); the more general case of correlated 
Gaussian coordinates can be treated by using a simple 
transformation of the coordinate system. As illustrated 
in Fig. 6, the elliptical area defined by the equation 



<xi v2tt 



\aj \a 2 J 



(3.2b) 



and 



-k 2 /2\ 



Py (P)=— Ue"^ 



<x 2 V 2tt 



><£<< 



(3.1b) 



where <xi and <x 2 are, respectively, the standard devia- 
tions of the x and y coordinates. Without loss of general- 



contains 100 (1 — e ) percent of the mobile positions, 
or about 39 % of the mobile positions when k = 1, 86 % 
when k = 2, and 99 % when k=3. The elliptical area 
containing nearly all the positions corresponds to the 
rectangular area shown in Fig. 1, so that the Gaussian- 
coordinate model can easily be related to the uniformly 
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1 i 

3a, 


39%v / s^~ 


^\ 66% 

1 ^ ry 


\V 


I Gj J2a, 30J 


v 


/ 99% 



Fig. 6. Elliptical areas associated Gaussian mobile coordinates. 
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distributed mobile model when it is convenient. For 
example, an ellipse just fitting inside the rectangle of 
Fig. 1 has the area \tsD\D 2 and contains \n = 78.54 % 
of the mobile positions for the rectangular, uniform dis- 
tribution. This same percentage for the Gaussian-coor- 
dinate model is contained in the elliptical area given by 
Eq. (3.2b) with k= 1.754, so that the two models are 
roughly equivalent when \ D\ ~ 1 .75 <xi and 
\ D 2 ~ 1.75 <r 2 , or Z)] ~ 3.5 <xi and D 2 ~ 3.5 <x 2 . 

Since a difference of independent Gaussian RVs with 
variances a and b is also a Gaussian RV whose variance 
is a + b, the differences in the coordinates of two mo- 
biles are Gaussian: 

Ax = Xi — Xj = G(0,2<x 2 ) and 

Ay = y;-y;=G(0,2cri) (3.3) 

where G(fJL, a 2 ) denotes a Gaussian RV with mean fju 
and variance <x 2 . The joint pdf of the differences is 
given by 



Pii,A y (a, P) ■ 



1 



4tTO"i(T2 



exp 



a 2 P 2 



. (3.4) 



The cumulative probability distribution function for the 
distance between two mobiles is formulated in terms of 
the squares of the Gaussian RVs A x and Ay as 

F,(y) = Pr{d(, < y } = Pr{\/(Ax) 2 + (Ay) 2 < y}. 

(3.5a) 



TT<Xi<X 2 
P 



2tht\(t 2 



2<X]<X 2 



/•ir/2 

/ dcf> exp{ — p\a + b cos 2cf>)} 
Jo 

J da exp{ —p 2 (a + Z? cos a)} (3.6b) 

Jq 



-ap* 



Ubp 2 ) 



(3.6c) 



in which we use the integral in Ref. [5], Sec. 9.6.16 to 
identify I (*)> the modified Bessel function of the first 
kind, and we define 



For convenience of notation and ease of comparison of 
the rectangular and Gaussian deployment models, we 
define the normalized variable £ = plD\ = p/Ka u 
where k = Dilat relates the dimensions of the rectangu- 
lar deployment area to the standard deviation of the 
Gaussian deployment distribution, and we denote the 
area shape parameter by £ = D\ID 2 = o-^a 2 to be consis- 
tent with the use of this symbol for the rectangular 
deployment area. Then the pdf of the link distance can 
be written 



• ^e-* 2 * 14 ^ loU 2 e(\ - £ 2 )/A P ^ (3.7a) 



Let us define the rectangular-to-polar change of vari- 
ables given by A x = dg cos and Ay = dg sin 0. The 
joint pdf of dij and 0, expressed in terms of the dummy 
variables p and <f>, is found to be 



P ^'^ ) = 4^ eXp 



cos 2 cft sin 2 cf> 

2 "^ 2~ 

. 0"! 0"2 



< <j> < 2tt, p > 0. 



(3.6a) 



The marginal pdf of d t j is found by integrating out the 
variable cf> in Eq. (3.6a). Noting that the joint density is 
the same in each of the four quadrants, we can write 

/•ir/2 

P</(p) = 4 / d<£p4*(p, 0) 



cos 2 sin 2 

2 ' 2 



TTO"i02. 



dcj) expj - ^ 



with the special case for £ = 1 (<T] = <r 2 ) given by 

(3.7b) 

Plots of the link distance pdf Eq. (3.7a) are shown in 
Fig. 7 for k = 3 (the length of the side of the rectangular 
deployment area is three times the standard deviation of 
the Gaussian deployment area in each direction) and 
£ = 1, 0.5, and 0.25. The similarity of these plots to the 
in Fig. 5 is strong; the similarity can be made even 
stronger by choosing a little smaller value than k = 
DJcTi = 3. Of course, the curves in Fig. 7 are smoother 
than those in Fig. 5 because the deployment area for the 
assumption of a Gaussian distribution of mobile loca- 
tions has no edges. 
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Now having the pdf of d tj , we can write the cdf Eq. 
(3.5a) for this RV as 

F d (y =K<r l ij)= / dpp^p) 



f Auu^ u2 ^ 2),% 1 {k 2 u\\ - £ 2 )l%\. (3.8a) 



■& 



For the special case of £= 1, Eq. (3.8a) becomes 
F d (y=K<T^) = Yj duue- Klu2/4 =J dv e" 



= 1 - e"* 2 ^ 4 . 



(3.8b) 



Plots of Eq. (3.8a) for k = 3, obtained by numerical 
integration, are shown in Fig. 8. 

3.2 Median and Mode for the Link Distance and 
Gaussian-Distributed Locations 

For £= 1, Eq. (3.7b) is easily differentiated to find the 
mode of the distribution and Eq. (3.8b) is easily solved 
for the median: 



£nod/ =1 = V2/K= 1.4142, (3.9a) 

£ me / =1 = 2Vln~2/K = 1.3863/k. (3.9b) 

From Table 2, the mode of the distribution for a random 
distribution of mobile locations in a rectangular area for 
£= 1 is 0.4786; the mode for the Gaussian distribution 
of mobile locations for £=1 matches it when 
k = 2.9549. From Table 1, the median of the distribution 
for a random distribution of mobile locations in a rectan- 
gular area for £= 1 is 0.5120; the median for the Gaus- 
sian distribution of mobile locations for £ = 1 matches it 
when k = 2.7076. 

4. Conclusions 

We have found the distributions for the distance be- 
tween randomly distributed mobiles for two different 
assumptions: (1) the mobile locations are uniformly dis- 
tributed in a rectangular area, and (2) the mobile loca- 
tions have a two-dimensional Gaussian distribution. The 
cdfs for both cases are very similar despite the fact that 
the first distribution has a finite boundary and the sec- 
ond does not. The implication of this finding is that for 
simulation or analysis of mobile communication sys- 
tems, the model used for the distribution of the mobile 
locations can be chosen for convenience. 
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Fig. 8. Plot of the link distance cdf for Gaussian- distributed mobile locations. 



5. Appendix A. Details of the Derivation 
of the Link Distance Distribution for a 
Rectangular Area 

The development here follows that in [2] but in more 
detail, correcting several typographical errors in that 
presentation. 



5.1 Evaluation for the Interval < y < D i9 or 

For this interval we use Eq. (2.6c); the upper limit of 
the first (outer) integral equals £ and t he upper limit of 
the second (inner) integral equals £V£ 2 — u 2 . Then, 

F,(y=£Z>i) = 4/ du(l-u) dv(l-v) 

•J Jo 

= 4f / dw(l -fw) / dv(l-v) 



= 4£f 2 / dwVl - w 2 - 4tf 3 J dw w VT^ 

*^0 *^0 



W 



= 2f/ dw(l-fw) 



l-^i-^vT^w 21 
= 2{ [ dw(l - fw)[2^Vl - w 2 - £ 2 £ 2 (1 - w 2 )] 



- 2£ 2 f /* dw(l - £w - w 2 + £w 3 ). (5.1a) 

From Ref. [3], Sec. 3.251.1 we have 

J dww^Ml -w A j =X B (a'V (5 * lb) 

where B(a, Z?) = r(a)r(^)/T(a + b) is the Beta func- 
tion. Applying Eq. (5.1b) to (5.1a) yields 

F,(y=^Z) 1 ) = 4^ 2 4B(i!)-4^ 3 4B(l,|) 

= # 2 [TT-5f(l+£)+3# 2 ]- (5.1C) 



5.2 Evaluation for the Interval Z)i < y < Z) 2 > or 

For this interval we use Eq. (2.6c); the upper limit of 
the first (outer) integral equals 1 and t he upper limit of 
the second (inner) integral equals £V£ 2 — u 2 . Then, 
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F d (y=^D,) = 4 I dw(l - u) I dv(l - v) = 1-2 / dw(l - u) I 1 - £V£ 2 - w 2 



Vi2Zp 



2 / dw(l - w) 



£2 - M 2 - £ 2 (f - M 2 )] 



i - \\ - i\/ir~u 

= 2 f du(l -u)[2£Vg- 

= 4f / duVg 2 - u 1 -AU duuV% 

J J 

-2C l du(f - £ 2 u - u 2 + u 3 ). 
Jo 

From Ref. [4], integral No. 157 we have 
1 



(5.2a) 



duVg 2 - u 2 -- . 



iVf^^+^sin-^l 



(5.2b) 



and in Ref. [4], integral No. 162 we have 



du uV? - u 2 = - ^ ( ¥ ~ u 2 



(5.2c) 



Substituting Eqs. (5.2b) and (5.2c) in Eq. (5.2a), we 
obtain 



F rf (y=££>,) = 2£ 



VF^T + fsin-'d 



U 



e - (e - 1) 



i\¥ - z) 



= i zvw^ve + 1) + 2tf 2 sin- i (i 



la*e 3 +6& 2 -a (s-w 



5.3 E valuation for the Interval D 2 ^ y ^ 
yJDl + Dl, or f" 1 < ^< VTTr 3 

For this third interval we use Eq. (2. 7c); the lo wer 
limit of the first (outer) integral equals V£ 2 — £ -2 and 
the lower li mit of the second (inner) integral equals 
£V£ 2 - w 2 . Then, 

F„(y=fZ) 1 )=l 
- 4 / d«(l - «) / dv(l - v) (5.3a) 






1 -2 / dw(l - u) 



1 - 2^V£ 2 - u 2 



= 1-2 



+ t 2 (e-u 2 )] 

/ dw(l -u) + 4U dwV^ 2 - u 2 



l ^e z F 1 



l ^e z F 1 



4f / du uVg 2 - u 2 



2£ : 



JVl 



dw(£ 2 - i; 2 U — U 2 + U 3 ) 



vprpi 



= i-(i-vf^t) 2 



+2 ^ vF^+^w^-r'V^r 
+^[(^ 2 -i) 3/2 -r 3 ] 



2^ 2 -r-i+i 



- ^ 2 v^ 2 -r 2 +^ 2 « 2 - r 2 )+kf 2 - r 2 ) 
-ae-c 2 ) 2 ] 

=i^V^T(2^ 2 + l) + 2^ 2 {sin- 1 (i)-cos- 1 (^)} 



'Vf^i 



* c^e^p- 



+ |V^ 2 - rWf 2 + 1) - K 2 (£ 4 + 2^ 2 - i) 
+ k" 2 "^ (53b) 

5.4 Special Case: D x =D 2 or f = 1 

When A = D 2 = D or £ = 1 , the interval from £ = 1 to 
£= £ _1 vanishes, and from Eqs. (5.1c) and (5.3b) the 
cumulative probability distribution for the distance be- 
tween any two mobiles becomes Eq. (5.4). A form of 
the result for this special case was published in [1]. 
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F d (y=fD) = 



fW "!£+*), 



£<0 
0<£<1 



ve^(2e+\)-{\e+u 2 -\) 

+ 2£ 2 [sirT 1 (l/£) - cos _1 (l/f), 1 < £< V2 

1, £^V2 



(5.4) 
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